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We study the dynamics of the the spin zero open charm and beauty mesons using QCD spectral sum rules 
(QSSR), where we observe the important role of the chiral condensate (ijjip) in the mass-splittings between 
the scalar-pseudoscalar mesons. Fixing the sum rule parameters for reproducing the well-known D(0~) and 
D s (0~) masses, we re-obtain the running charm quark mass: fh c (m c ) = 1.131q'q4 GeV, which confirms 
our recent estimate from this channel P^. Therefore, using sum rules with no-free parameters, we deduce 
£f~l • M I ) 3 ( +) — (2297 ± 113) MeV, which is consistent with the observed D s (2317) meson, while a small SU(3) 

J> ' breaking of about 25 MeV for the D s (0 + ) — D(0 + ) mass-difference has been obtained. We extend our analysis 

to the i3-system and find M B ( +j — -Ms(o-) — (422 ± 196) MeV confirming our old result from moment sum 
rules 0. Assuming an approximate (heavy and light) flavour and spin symmetries of the mass-splittings 
L^J . as indicated by our results, we also deduce M^.^+j ~ (2440 ± 113) MeV in agreement with the observed 

L> sJ (2457). We also get: f D (o+) - (217 ± 25) MeV much bigger than ^=130.6 MeV, while the size of the 
£f) . SU(3) breaking ratio /_d„(o+)//d(o+) — 0-93 ± 0.02 is opposite to the one of the 0~ channel of about 1.14. 
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1 Introduction 



Ph! The recent BABAR, BELLE and CLEO observations of two new states L> sJ (2317) and L> sJ (2457) 3 in the 

D s ir, -D s 7 and D s ttj final states have stimulated a renewed interest in the spectroscopy of open charm states 
which one can notice from different recent theoretical attempts to identify their nature [3]. In this paper, 
we will try to provide the answer to this question from QCD spectral sum rules a la Shifman-Vainshtein- 
Zakharov [4. In fact, we have already addressed a similar question in the past [2], where we have predicted 
using QSSR the mass splitting of the + — 0~ and l~ — 1 + bu mesons using double ratio of moments sum 
rules based on an expansion in the inverse of the b quark mass. We found that the value of the mass- 
splittings between the chiral multiplets were about the same and approximately independent of the spin of 
these mesons signaling an heavy quark- type approximate symmetry: 

Mb(o+) ~ Mb (0 -) » M B , {1+) - Mb.(i-) ps (417 ± 212) MeV . (1) 

The effect and errors on the mass-splittings are mainly due to the chiral condensate (i/jip) an d to the value 
of the b quark mass. In this paper, we shall use an analogous approach to the open charm states. However, 
the same method in terms of the 1 /m c expansion and some other nonrelativistic sum rules will be dangerous 
here due to the relatively light value of the charm quark mass. Instead, we shall work with relativistic 
exponential sum rules used successfully in the light quark channels for predicting the meson masses and 
QCD parameters [H] and in the D and B channels for predicting the (famous) decay constants Jd.b P3IE1IE] 
and the charm and bottom quark masses El El Ej ■ 

2 The QCD spectral sum rules 

We shall work here with the (pseudo)scalar two-point correlators: 

^p/s(q 2 ) = *Jd 4 x e**" (0\TJ P/S (x)4 /S (p)\0), (2) 
built from the (pseudo)scalar and (axial)-vector heavy-light quark currents: 



Jp/s{x) = (m Q ± m q )Q(i-f 5 )q, J^ /A = Q^{^)q 



(3) 



If we fix Q = c and q = s, the corresponding mesons have the quantum numbers of the D s (0~), D s (0 + ) 
mesons, vtlq and m s are the running quark masses. In the (pseudo)scalar channels, the relevant sum rules 
for our problem are the Laplace transform sum rules: 



dt e 



lm.ippi s {t), and 



dr 



log Lf /S {r) 



(4) 



where t< is the hadronic threshold, and H denotes the corresponding meson. The latter sum rule, or its 
slight modification, is useful, as it is equal to the resonance mass squared, in the simple duality ansatz 
parametrization of the spectral function: 



1 



lmt/jp (t) ~ fz)Mf)5(t — Mp) + "QCD continuum" 6 i c ), 



(5) 



where the "QCD continuum comes from the discontinuity of the QCD diagrams, which is expected to give 
a good smearing of the different radial excitations . The decay constant fp> is analogous to f„ = 130.6 
MeV; t c is the QCD continuum threshold, which is, like the sum rule variable t, an (a priori) arbitrary 
parameter. In this paper, we shall impose the r- and ^-stability criteria for extracting our optimal results. 
The corresponding t c value also agrees with the FESR duality constraints HI and very roughly indicates 
the position of the next radial excitations. However, in order to have a conservative result, we take a largest 
range of t c from the beginning of r- to the one of ^-stabilities. 

The QCD expression of the correlator is well-known to two-loop accuracy (see e.g. [S] and the explicit 
expressions given in in terms of the perturbative pole mass Mq, and including the non-perturbativc 

condensates of dimensions less than or equal to six 2 . For a pedagocial presentation, we write the sum 
rule in the chiral limit (m s = 0), where the expression is more compact. In this way, one can understand 
qualitatively the source of the mass splittings. The sum rule reads to order a s : 
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The different terms are 
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where we have used the contribution of the gluon condensate given in which is IR finite when letting 
4 . The previous sum rules can be expressed in terms of the running mass ffiQ(v) through the 



perturbative two-loop relation 

Mq = m Q {p 2 ) 




(9) 



where Mq is the pole mass. Throughout this paper we shall use the values of the QCD parameters given in 
Table 1. 



1 At the optimization scale, its effect is negligible, such that a more involved parametrization is not necessary. 
2 We shall include the negligible contribution from the dimension six four-quark condensates, while we shall neglect an 
eventual tachyonic gluon mass correction term found to be negligible in some other channels J10| , 

3 Notice a missprint in the expression given in which does not affect the result obtained there. 

4 The numerical change is negligible compared with the original expression obtained in [121 . 



Table 1: QCD input parameters used in the analysis. 



Parameters References 



A 4 = (325 ± 43) MeV 

A 5 = (225 ± 30) MeV 

fh b (m b ) = (4.24 ±0.06) GeV 00 Ed 

m s (2 GeV) = (111 ± 22) MeV El EE] 

(dd) 1 / 3 (2 GeV)=-(243± 14) MeV 001111 

(ss)/{dd) = 0.8 ±0.1 0EFJ 

{a s G 2 } = (0.07 ±0.01) GeV 4 EH] 

Af 2 = (0.8 ± 0.1) GeV 2 00 

a s (il}iji) 2 = (5.8 ± 2.4) x 10~ 4 GeV 6 |S1 CU EDJ 



We have used for the mixed condensate the parametrization: 

g(da^G^d) = M 2 (dd), (10) 

and deduced the value of the QCD scale A from the value of a s {M z ) = (0.1184 ± 0.031) H21 We have 
taken the mean value of m s from recent reviews [51 181 IH)]. 



3 Calibration of the sum rule from the D(0 ) and D s (0 ) masses 
and re-estimate of m c (m c ) 




Figure 1: r in GeV~ 2 -dependence of the a) Md(0") in GeV for m c (m c ) — 1.11 GeV and b) M Da (o-) in 
GeV for fh c (m c ) = 1.15 GeV at a given value of t c = 7.5 GeV 2 . The dashed line is the result including the 
leading (ifnp) contribution. The full line is the one including non-perturbative effects up to dimension-six. 



This analysis has been already done in previous papers to order a s and a 2 and has served to fix the running 
charm quark mass. We repeat this analysis here to order a s for a pedagogical purpose. We show in Fig. 
la), the r-dependence of the D(0~) and in Fig lb) the one of the D s (0~) masses for a given value of t c , 
which is the central value of the range: 

t c = (7.5 ± 1.5) GeV 2 , (11) 

where the lowest value corresponds to the beginning of r-stablity and the highest one to the beginning of 
t c stability obtained by ^ El m the analysis of fn and fo s ■ This range of t c -values covers the different 
choices of t c used in the sum rule literature. As mentioned previously, the one of the beginning of testability 
coincides, in general, with the value obtained from FESR local duality constraints 00. Using the input 



values of QCD parameters in Table 1, the best fits of D(0~) (resp. D s (0~)) masses for a given value of 
t c = 7.5 GeV 2 correspond to a value of m c (m c ) of 1.11 (resp. 1.15) GeV. Taking the mean value as an 
estimate, one can deduce: 

fh c {m 2 c ) = (1.13±8;8£ ± 0.02 ± 0.02 ± 0.02) GeV , (12) 

where the errors come respectively from t c , (i/V), A and the mean value of m c required from fitting the 
D(0~) and D s (0~) masses. This value is perfectly consistent with the one obtained in P3 E] obtained to 
the same order and to order a 2 , indicating that, though the a 2 corrections are both large in the two-point 
function and m c [23], it does not affect much the final result from the sum rule analysis. In fact, higher 
corrections tend mainly to shift the position of the stability regions but affect slightly the output value of 
m c . This value of m c is in the range of the current average value (1.23 ± 0.05) GeV reviewed in [51 151 |2"T]. 
However, it does not favour higher values of m c allowed in some other channels and by some non relativistic 
sum rules and approaches. However, these non relativistic approaches might be quite inaccurate due to the 
relative smallness of the charm quark mass. Higher values of m c would lead to an overestimate of the D(0~) 
and D 3 (0~) masses. In the following analysis, we shall use the central value fh c (m c ) = 1.11 (resp. 1.15) 
GeV for the non-strange (resp. strange) meson channels. 

4 The + — meson mass-splittings 
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Figure 2: Similar to Fig. 1 but r-behaviour of Mr> s (o+) for given values of t c — 7.5 GeV 2 and m c {m c ) = 1.15 
GeV. 



• We study in Fig. 2), the r-dependence of the D s (0 + ) mass at the values of t c and m c obtained 
previously. In this way, we obtain: 

M Ds{Q+) ~ (2297±|8^q±11±2±11) MeV =S> M Db[q+) - M^jo-j = (328± 113) MeV , (13) 

where the errors come respectively from t c , m c , (tpip), m s , and A. We have used the experimental value 
of M Ds (Q-y The reduction of the theoretical error needs precise values of the continuum threshold 5 
and of the charm quark mass which are not within the present reach of the estimate of these quantities. 
Further discoveries of the continuum states will reduce the present error in the splitting. One should 
also notice that in the ratio of sum rules with which we are working, we expect that perturbative 
radiative corrections are minimized though individually large in the expression of the correlator and 
of the quark mass. 

• The value of the mass-splittings obtained previously is comparable with the one of the B(0 + )-B(0~) 
given in Eq. ft}, and suggests an approximate heavy-flavour symmetry of this observable. 

• We also derive the result in the limit of SU(3)f symmetry where the strange quark mass is put to 
zero, and where the (tptp) condensate is chirally symmetric ((ss) = (dd)). In this case, one can predict 
an approximate degenerate mass within the errors: 

M Ds{0+) - M D{0+) ~ 25 MeV , (14) 

5 The range of t c values 6-9 GeV 2 obtained previously for the D(0~ ) mesons coincides a posteriori with the corresponding 
range for the D(0 + ) meson if one assumes that the splitting between the radial excitations is the same as the one between the 
ground states, i.e about 300 MeV. We have cheked during the analysis that this effect is unimportant and is inside the large 
error induced by the range of t c used. 



which indicates that the mass-splitting between the strange and non-strange + open charm mesons 
is almost not affected by SU(3) breakings, contrary to the case of the 0~ mesons with a splitting of 
about 100 MeV. 

• We extend the analysis to the case of the B(0 + ) meson. Here, it is more informative to predict the 
ratio of the + over the 0~ masses as the prediction on the absolute values though presenting stability 
in r tend to overestimate the value of Mb ■ We obtain: 

( ' ~ (1.08 ±0.03 ±0.03 ±0.02 ±0.02) M B{0+) - M B{0 -) ~ (422 ± 196) MeV , (15) 

where the errors come respectively from t c taken in the range 43 — 60 GeV 2 , m&, (^ip), and r. It agrees 
with the result in Eq. QJ obtained from moment sum rules We have used the value of fhb{mb) 
given in Table 1. 

5 The 1 + — 1 meson mass-splitting 

Our previous results in Eqs. QJ, ifH5|) to (|l"5|) suggest that the mass-splittings are approximately (heavy 
and light) flavour and spin independent. Therefore, one can write to a good approximation the empirical 
relation: 

Md s (o+)— Mn s (o-) ~ Md(o+)- m d(o-) ~ ^b(o+) _ -^b(o-) ~ Mb(i+)~Mb(i-) ~ ^_d*(i+)- a ^d*(i-) ■ (16) 
Using the most precise number given in Eq. I)13|l . one can deduce: 

M D*(i+) ~ M d;(i-) - ( 328 ± 113 ) MeV => M Dt{x+) ~ (2440 ± 113) MeV . (17) 

This result is consistent with the 1 + assignement of the cs meson D s j (2457) discovered recently [3]. In a 
future work, we plan to study in details this spin one channel using QSSR. 

6 The D(0 + ) and D s (0 + ) decay constants 




Figure 3: Similar to Fig. 1 but r-behaviour of /d(0 + ) for given values of t c = 7.5 GeV 2 and m c (m c ) = 1.11 
GeV. 



For completing our analysis, we estimate the decay constant Jd(o+) analogue to / ff = 130.6 MeV. We show 
the behaviour of /r>(o+) versus t, where a goood stablity is obtained. Adopting the range of £ c -values 
obtained previously and using fh c (m c ) = I.II^q'q^ GeV required for a best fit of the non strange D(0 + ) 
meson mass, we deduce to two-loop accuracy: 

^(0+) = (217±*+± 5 19 ± 10 ± 10) MeV , (18) 

where the errors come respectively from the values of t c , m c , ('tpip) and A. We have fixed M £ j( 0+ - ) to be about 
2272 MeV from our previous fit. It is informative to compare this result with the one of fa = (205 ± 20) 
MeV, where the main difference can be attributed by the sign flip of the quark condensate contribution in 



the QCD expression of the corresponding correlators. A numerical study of the SU(3) breaking effect leads 
to: 

_ /d,(o+) 



0.93 ±0.02 



D(0+) 



(19) 



which is reverse to the analogous ratio in the pseudoscalar channel JdJ Jd — 1.14 ± 0.04 given semi- 
analytically in [24 j . In order to understand this result, we give a semi-analytic parametrization of this SU(3) 
breaking ratio. Keeping the leading term in m s and (i^ip}, one obtains: 



1 - 



1 - 7.5(ss - dd) 



1/2 (M 



l D s (0+) 



D(0+ 



0.9 , 



(20) 



where the main effect comes from the negative sign of the m s contribution in the overall normalization of 
the scalar current, while the meson mass ratio does not compensate this effect because of the almost equal 
mass of D s (0 + ) and -D(0 + ) obtained in previous analysis. This feature is opposite to the case of fo(o-)- 



7 Summary and conclusions 

Due to the experimental recent discovery of the D s j(2317) and D* 7 (2457), we have analyzed using QSSR 
the dynamics of the and ^ open charm and beauty meson channels: 

• We have re-estimated the running charm quark mass from the D and D s mesons. The result in Eq. 
(|12fl confirms earlier results obtained to two- and three-loop accuracies 

• We have studied the mass-splittings of the + -0~ in the D systems using QSSR. Our result in the (0 + ) 
channel given in Eq. (|13|) agrees with the recent experimental findings of the D s j(2317) suggesting 
that this state is a good candidate for being a cs + meson. 

• We also found, in Eq. 1|14|) . that the SU(i) breaking responsible of the mass-splitting between the 
D s (0 + ) and D(0~) is small of about 25 MeV contrary to the case of the pseudoscalar D s -D mesons of 
about 100 MeV. 

• We have extended our analysis to the B-system. Our results in Eqs. JT3J, i|15|) and Q suggest an 
approximate (light and heavy) flavour and spin symmetries of the meson mass-splittings. We use this 
result to get the mass of the cs D*(l + ) meson in Eq. (|17|) . which is in (surprising) good agreement 
with the observed £^(2457). 

• Using QSSR, we have also determined the decay constants of the + mesons and compare them with 
the ones of the 0~ states. The result in Eq. i|18[l . which is similar to the pseudoscalar decay constant 
fr> — 205 MeV, suggests a huge violation of the heavy quark symmetry l/VA/jj scaling law. Finally, 
our results in Eqs. (|19|l and l|20|) indicate that the SU(3) breaking act in an opposite way compared 
to the case of the 0~ channels. 

Experimental or/and lattice measurements of the previous predictions are useful for testing the validity of 
the results obtained to two-loop accuracy in this paper from QCD spectral sum rules. However, a complete 
confirmation of the nature of these new states needs a detail study of their production and decays. We plan 
to come back to these points in a future work. 

During the editorial preparation of this paper, there appears in the literature a paper using sum rules 
method to the same channel [25]. Instead of our result in Eq. (13): M£> s ( +) = (2297 ± 113) MeV, the 
authors obtain M Ds (o + ) = (2480 ±30) MeV which is higher than the BABAR result D S (2317) by about 160 
MeV. Considering this deviation as significant, the authors conclude that the experimental candidate cannot 
be a cs state contrary to the conclusion reached in the present paper. However, by scrutinizing the analysis 
of Ref. [25], we find that the true errors of the analysis have been underestimated: 

• As one can see from their figures, the quoted error of 30 MeV only takes into account the one due to the 
choice of the QCD contiuum threshold taken in a smaller range 8.1-9.3 GeV 2 , than the more conservative 
value 6-9 GeV 2 used in the present paper, which induces a larger error of 80-98 MeV. 

• The high central value obtained in [25] is related to a higher choice of the (ill-defined) charm quark pole 
mass of 1.46 GeV compared to the value 1.3 GeV which would have been deduced from its relation with the 
running charm quark mass 1.15 GeV in Eq. (9) used in this paper to reproduce the D s (0~) mass. As the 
analysis is very sensitive to m c , its induced error should be included in the final number. With m c in Eq. 
(12), this uncertainty is about 70 MeV, and might be bigger if one considers all range of m c given in the 



literature. In this paper, we have used Mp andMo s as alternative channels for extracting m c . 
• Taking properly the different sources of errors including the running of condensate and the effects of some 
other anomalous dimensions, one then leads to a consistency within the errors of both theoretical estimates 
with the present data. 
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